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Polarization based sensing with tilted fiber Bragg grating (TFBG) sensors is analysed theoretically 
by two alternative approaches. The first method is based on tracking the grating transmission for 
two orthogonal states of linear polarized light that are extracted from the measured Jones matrix or 
Stokes vectors of the TFBG transmission spectra. The second method is based on the measurements 
along the system principle axes and polarization dependent loss (PDF) parameter, also calculated 
from measured data. It is shown that the frequent crossing of the Jones matrix eigenvalues as a 
function of wavelength leads to a non-physical interchange of the calculated principal axes; a method 
to remove this unwanted mathematical artefact and to restore the order of the system eigenvalues 
and the corresponding principal axes is provided. A comparison of the two approaches reveals that 
the PDF method provides a smaller standard deviation and therefore lower limit of detection in 
refractometric sensing. Furthermore, the polarization analysis of the measured spectra allows for 
the identihcation of the principal states of polarization of the sensor system and consequentially for 
the calculation of the transmission spectrum for any incident polarization state. The stability of the 
orientation of the system principal axes is also investigated as a function of wavelength. 


I. INTRODUCTION 

Polarization-based sensing is crucial for various types 
of optical sensors, in particular for stress analysis, 
plasmon-mediated sensing, and sensing of anisotropic 
media or other forms of perturbations [l|-Q. In the par¬ 
ticular case of waveguide-type sensors (including opti¬ 
cal fibers), it is possible to interrogate the device op¬ 
tical properties with polarized light but in general this 
require very careful alignment and control of the input 
polarization (which is relatively easy in waveguides but 
not so much in conventional fibers). Here we describe a 
data analysis method that provides a full characterization 
of the polarization-resolved optical transmission state of 
any sensor system based on measurements of the wave¬ 
length dependence of the Jones matrix with an Optical 
Vector Analyser (OVA). The OVA measures the trans¬ 
mission amplitude and phase (or reflection) of an optical 
system as a function of polarization and wavelength and 
provides the Jones matrix elements for each wavelength. 
We use a standard tilted fiber Bragg grating (TFBG) 
refractometric sensor as a test system for demonstrat¬ 
ing the new approach. The results further confirm that 
the polarization dependent loss parameter (PDF) of the 
system provides a smaller standard deviation (and hence 
better limit of detection) than the transmission spectra 
extracted along the principal axes, as initially observed 
empirically in [ 3 . 

It is known that TFBGs have strong polarization- 
dependent properties [6|-[8| due to the tilt of the grating 
planes which breaks the cylindrical symmetry of the fiber 
and strongly impacts the magnitude of the coupling co- 
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efficients between the incident HEn core mode and the 
cladding modes of various polarizations (TE, TM, HE, 
EH modes) excited by the grating [9| , as is schematically 
shown in Fig. [TJ 





FIG. I. The schematic representation of TFBG grating and 
the incident linearly polarized core mode (the Ep component 
is shown). The core mode is rotated about the optical device 
axis by some angle 0. The grating is tilted by 0 angle about 
the X axis. 

In particular, for a grating tilted by an angle 0 in the 
reference y — z plane and a linearly polarized input core 
mode whose polarization is rotated by an arbitrary an¬ 
gle (j) with the reference x axis (a most general occurrence 
in optical fibers that are not polarization preserving), 
the coupling to individual cladding modes will depend 
strongly on both 0 and (f. It has been further demon¬ 
strated recently that high order cladding modes of the 
EH and TM families have radially polarized evanescent 
fields around the fiber cladding while HE/TE modes 
are azimuthally polarized. This is important because in 
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many sensing applications (or other forms of interactions 
between cladding guided light and coatings or external 
media, such as nonlinear switching 0 , or four wave 
mixing applications), the state of polarization of the 
evanescent field can be a determinant factor. This most 
easily seen for metallic coatings whose boundaries condi¬ 
tions depend strongly on the polarization of the incident 
electric field relative to the plane of the boundary [El, 0 . 
It was also demonstrated in [H that only radially po¬ 
larized cladding modes are excited by TFBGs when the 
input core mode is linearly polarized along the tilt plane 
{i.e. along ^), while the orthogonal input polarization 
(along x) generates only azimuthally polarized cladding 
modes. For convenience and for easy comparison with 
conventional optics, in the remainder of this paper we 
use ”TE” to refer to the HE and TE fiber mode families 
(azimuthally polarized, i.e. polarized along the plane of 
the fiber cladding boundary) and ” TM” for the EH and 
TM fiber mode families (radially polarized, i.e. perpen¬ 
dicularly to the fiber cladding surface). This is further 
justified by the fact that on the scale of the wavelength 
(1.5 /im, the curvature of the cladding boundary (diam¬ 
eter of 125 iim) is negligible and the boundary can be 
viewed as approximately flat surface. The effect of the 
input mode polarization can be observed by measuring 
optical transmission spectra with a polarizer inserted be¬ 
tween the light source and the grating [T2|. A typical 
series of spectra measured at various angles of linearly 
polarized incident light is shown in Fig.[2]for a 1 cm long 
TFBG immersed in water. 

It is clear from Fig. [2] that there are two orthogonal 
states of transmission, oriented at 90^ degrees relative 
to each other, where the spectrum has well-defined ex¬ 
trema. We can therefore find the polarizer angles for 
which TE and TM modes are excited a posteriori with¬ 
out prior knowledge of the actual orientation of the tilt 
plane. However this is often impractical, especially if 
the device under test has to be moved between measure¬ 
ments, thereby potentially altering the light polarization 
orientation between the polarizer and the TFBG. A much 
more practical and rapid method to interrogate such de¬ 
vices is to use an OVA to get the Jones matrix elements 
at each wavelength and to use the information to calcu¬ 
late the full set of polarization properties of the system, 
including the transmission along the chosen geometrical 
axes {Ix and ly in Fig. [2]) as well as the transmission along 
the system principal axes and PDF spectrum. How this 
is done is presented in the next section. 


II. MEASUREMENTS ALONG PRINCIPAL 
AXES OF AN OPTICAL SYSTEM 

In this section we review the PDL technique and pro¬ 
vide an approach to study the system transmission spec¬ 
trum along its principal axes. 

A linear optical system can be represented in terms of 
the 2x2 Jones matrix [J] connecting the incident and 


1(90°) = ly 
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FIG. 2. A typical TFBG transmission spectrum for linearly 
polarized light, and series of spectra obtained by rotating a 
linear polarizer about the optical axis are shown as the density 
plot. 

transmitted electric field vectors [T3| : 

Eout = [J]Ein ( 1 ) 

In practice, the transmission spectrum of a device un¬ 
der test is usually characterized by two major param¬ 
eters: a polarization-independent parameter called the 
insertion loss /(A) and the polarization-dependent loss 
parameter PiJL(A), both measured as a function of the 
free space wavelength A. To extract these parameters 
from the Jones matrix, a hermitian matrix [H] has to be 
constructed first. 

[H] = [J]t[J] (2) 

Next, carrying out the eigen decomposition of the [H] 
matrix 

[H] = [U][A][ur (3) 

allows us to find a diagonal matrix [A] containing the 
eigenvalues pi and p 2 of the [H] matrix. 

Alternatively, the singular value dec omp osition (SVD) 
of the Jones matrix can be computed [l^ : 

[j] = [umivf (4) 

where the diagonal matrix [S] contains two singular val¬ 
ues (Ji and (72 such that pi = and p 2 = cr 2 ^ due to the 
fact that [H] = 
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The matrix [U] = resulting from the SVD con¬ 

tains two orthogonal vectors corresponding to pi and p 2 
eigenvalues, and these eigenvectors correspond to the de¬ 
vice response when the electric field vector E is aligned 
with the system principal axes, i.e. E \\ uj^ of the system. 
The eigenvalues pk = can be interpreted as the ob¬ 
servable transmission loss for the corresponding incident 
polarizations. The principal axes fli, U 2 of an arbitrary 
fiber device are showed schematically in Fig. IH The an¬ 
gle (j) corresponds to the rotation between the arbitrarily 
selected frame of reference (x, y) and the principal axes 
defined by fli, 



FIG. 3. TFBG with physical x, y axes and principal system 
axes ui, U 2 measured at some wavelength A. 

Having the singular values of the Jones matrix [J] or 
eigenvalues of [H] = [J]^[J] matrix allow for the compu¬ 
tation of the insertion loss, which is the average loss over 
all possible polarizations, as follows HI. 

lout = {Eout\Eout) = {Ein\[J]^[J]\Ein) = 

= {E'inlME'in) = ( 5 ) 

Which can be re-formulated in the more frequently used 
dB scale as: 

lout plEp2 

— = 101ogio—-—, (6) 

J-in ^ 

here, the scalar product is denoted in accordance with 
Dirac’s braket notation. Similarly, the polarization de¬ 
pendent loss, which is the magnitude of the difference be¬ 
tween the maximum and minimum device transmission 
over all possible input polarization actually corresponds 
to the difference in the loss measured along the system 
principal axes ui and U 2 . It is defined as follows: 

Pi?L = 10|logio-|. (7) 

P2 

Alternatively we can introduce a polarization parame¬ 
ter with a linear scale Q: 


In addition to the eigenvalues of the system trans¬ 
mission matrix (plotted as a function of wavelength in 
Fig. 0]), we can also find the angle of rotation (0) of the 


principal axes about the device optical axis {z) and plot 
it as a function of wavelength (bottom frame of Fig. 0]). 
From the mathematical point of view, both eigenvalues 
are roots of quadric polynomial and can be order ar¬ 
bitrary, usually in descending order: this is why pi is 
always larger than p 2 in Fig. 0^, hence the eigenstates 
are interchanged when they cross. Because of this ef¬ 
fect, the corresponding eigenvectors are interchanged as 
well, therefore the angle 0 experiences 7r/2 shifts at the 
crossing points, as shown on the bottom panel of Fig. 0^. 
Therefore, a strategy for restoring the individual trans¬ 
mission spectra along the principal axes is to interchange 
the eigenvalues and principal axes every time the 0 = 7r/2 
jumps are detected. 

The result of such reordering is shown Fig. 0)3. Now 
the transmission spectrum corresponds to linearly polar¬ 
ized light with its electric field vector aligned with each 
principal axis of the system and the rotation angle jumps 
are eliminated. 

In general, the principal axes of an optical system are 
not necessarily fixed with respect to a reference frame 
but may depend on the wavelength, as shown in Fig. [5] 
for a TFBG sensor. Indeed, by changing the optical fre¬ 
quency it is expected that an optical system would op¬ 
erate differently. The global behavior of the principal 
axes, along the whole operational range of the TFBG 
sensor, is shown in FiglSb, and the corresponding inser¬ 
tion loss in Figl5^. To remove the noise only the points 
at which the difference between eigenvalues is noticeable 
{\p 2 — pi\> 0.05) are plotted. The noise arises from the 
fact, that at the points of eigenvalues crossing the trans¬ 
mission matrix become degenerate, has the eigenvectors 
are not well defined, and oscillate rapidly with respect to 
the geometrical axis, as it can be seen from Fig. [Sb at 
the points of crossing. 

Figure [SJ) shows that the system possesses significant 
polarization asymmetry, or birefringence, in the wave¬ 
length range A G [1545 — 1575]. Although the princi¬ 
pal axes are globally stable, near 75*^ degrees relative 
to the reference frame of the LUNA interrogation sys¬ 
tem, locally they experience small oscillations, of about 
8^ degrees about the optical axis. These oscillations 
are related to the resonances observed in the spectrum 
(Fig. [5^), and reflect a wavelength dependent birefrin¬ 
gence. 

We also note the observed alternation between the 
peaks (Fig. [Sb), which can be explained by the fact that 
the alternating peaks have different azimuthal symme¬ 
tries and polarizations [9[. 


III. EXTRACTING TRANSMISSION SPECTRA 
FOR LIGHT POLARIZED IN AN ARBITRARY 
DIRECTION FROM THE JONES MATRIX OR 
THE STOKES VECTOR DATA 

In this section we describe an approach allowing to 
extract the transmission loss spectra for light that is lin- 
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FIG. 4. Eigenvalues pi and p 2 and the angle 0 between the geometrical axes of the system x, y and the coordinate system 
defined by the principal axes ui, U 2 , before a) and after b) the eigenvalues were reordered. (The data were obtained by means 
of the OVA 5000 Luna Technologies.) 
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FIG. 5. Rotation of the principal axes of TFBG device as a 
function of optical wavelength. (The data was obtained by 
means of the OVA 5000 Luna Technologies.) 


A. The Jones matrix analysis 

The Jones matrix of an optical device under test con¬ 
nected to a linear polarizer is given by: 


[Jd ev-\-pol] — [Jdev][Jpol]' 


(9) 


Here [Jdev] is Jones matrix of a device measured by an 
optical vector analyzer and [Jpoi] is the Jones matrix of 
a linear polarizer, oriented at an arbitrary angle theta 
relative to the system frame of reference given by: 


[jpoim = 


cos^(6>) cos{0) sin(6>) 
sin(6>) cos(^) sin^(^) ’ 


( 10 ) 


Thus, on optical device with a known Jones matrix 
[Jdev]^ connected to a linear polarizer rotated by 0 angle, 
is described by the following Jones matrix: 


[Jdev-\-pol{^)] — [Jde 


cos^(6>) cos{0) sm{0) 
sm{0) cos{0) 


( 11 ) 

And finally the average transmission spectrum /(^, A), 
for a given angle 0 of the linear polarizer, can be com¬ 
puted as 


early polarized in any direction, which may be a useful 
predictive tool for device performance under real life con¬ 
ditions without the need for an external polarizer. 


I{0, A) = 10 logio 


pi{ 0 , A) + p 2 { 0 , A) 
2 


( 12 ) 


where pi{0,X) and p 2 { 0 ,X) are the eigenvalues of 
H{0^XXj = \Jdev-\-poi{P 1 \Jdev-\-poi{P matrix. 
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B. The Stokes vector analysis 

As a supplementary information, we provide here an 
alternative approach to calculating the response of a de¬ 
vice to arbitrarily oriented linearly polarized light us¬ 
ing the Stokes vector of the system. The Stokes vector 
can be measured with a polarization controller or with a 
dedicated instrument, such as the JDS Uniphase SWS- 
OMNI-2 system. 

A beam of li ght can be completely described by four 
parameters [l5|, [l^ , represented in the form of the Stokes 
vector: 


following expressions: 

h{e, A) = i(5o(A) + Si{\) cos(20) + 52(A) sin(20)), 

= h5o(A) - 5i(A)cos(20) - 52(A) sin(20)). 

(18) 

Hence, if the Stokes vector is know in one coordinate 
system, the transmission spectrum of linearly polarized 
light with the electric field E aligned along an arbitrary 
angle 6>, can be recovered. 


S = 


( /(O^) + /(90^) \ 
/(O^) - /(90^) 
/(45^) - 1(135^) 

V ^RHS — IlHS / 


(13) 


IV. COMPARISON BETWEEN 
MEASUREMENTS ALONG THE 
GEOMETRICAL AND PRINCIPAL AXES 


here I{0) is the measured power transmission coefficient 
of the light polarized in the direction defined by the angle 
0 in the plane perpendicular to the direction of light prop¬ 
agation, and Irhs, Ilhc are the coefficients for right- 
and left-handed circular polarized light, respectively. 

Since the Stokes parameters are dependent upon the 
choice of axes, they can be transformed into a different 
coordinate system with a rotation matrix [R ]. 


[Rm = 


cos{0) 

sm{0) 


— sin(^) 
cos(<9) ’ 


(14) 


Considering that a second coordinate system is ob¬ 
tained by rotating the original coordinate system about 
the direction of light propagation on the angle 6>, we can 
write [13 


S'{0) = [R{0)]S, 


(15) 


or 


/S'o\ 

^2 

\sU 


TO 0 0" 

0 cos(2^) sin(2^) 0 

0 -sin(2<9) cos(2(9) 0 

_0 0 0 1 _ 



/ \ 

Si cos(26>) + S 2 sm{20) | 
—Si sm{20) + S 2 cos{20) I ’ 

\ Ss J 


(16) 


Now considering that 

1 ( 0 ") + 1 ( 90 ") = 

/(O^) - /(90^) = Si cos(2<9) + S 2 sin(2(9), (17) 


we conclude that the transmission loss spectra along the 
two orthogonal axes, rotated by the angle 0 with respect 
to the original system of measurements, are given by the 


The last results show that the transmission spectra of 
light linearly polarized along any direction in the trans¬ 
verse plane of the optical axis can be computed from the 
Jones matrix or from the Stokes vector, in addition the 
possibility of measuring it directly by rotating a polar¬ 
izer, as shown in Fig. [2l In particular, the transmission 
spectra of light linearly polarized along the system prin¬ 
cipal axes can be extracted from the Jones matrix and 
compared to the direct measurement with a polarizer. 
Figure [6] shows a comparison between the direct mea¬ 
surement of one resonance, obtained by aligning the po¬ 
larizer along the average direction of the principal axes 
(for a given spectrum), and the eigenstate spectra ob¬ 
tained from the Jones matrix. The small difference in the 
two results comes from the aforementioned wavelength 
dependent oscillation of the principal axes (which can¬ 
not be compensated for in the direct measurement), and 
to a possible slight change in polarization state between 
polarizer and the grating (since a non-polarization main¬ 
taining fiber is used). The following section will demon¬ 
strate that in spite of this small inaccuracy, the parame¬ 
ters extracted from the Jones matrix data provide excel¬ 
lent spectral sensitivity results for refractometric sensing. 


V. POLARIZATION-BASED DETECTION OF 
SMALL REFRACTIVE INDEX CHANGES WITH 
TFBG SENSORS 

In this section we investigate which polarization-based 
measurement techniques provide the best signal-to-noise 
ratio when TFBG sensors are used to detect small refrac¬ 
tive index changes, and use the special case of a TFBG 
coated with gold nanorods (as described in 0 )- This 
choice is made because the polarization dependence of 
waveguide-type sensors is much enhanced when metal 
coatings are used. As indicated earlier, the interaction 
of guided waves with metal interfaces depends strongly 
on whether the electric fields of the waves are tangential 
or normal to the metal boundary. It was further men¬ 
tioned that when the core-guided input light of a TFBG 
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FIG. 6. Transmission losses along TFBG system principal 
axes and its geometrical axes for perfectly aligned coordinate 
systems 0 = 0^. 


is linearly polarized along the principal axes, the elec¬ 
tric fields of high order cladding modes are either tan¬ 
gential or radial (hence normal) to the cladding bound¬ 
ary. To be precise, ^-polarized input light (correspond¬ 
ing to light polarized parallel to the plane of incidence 
on the tilted grating fringes, as shown in Fig. [U i.e. P- 
polarized) couples to radially polarized cladding modes 
while x-polarized light (perpendicular to the tilt plane, 
or S-polarized) couples to azimuthally polarized cladding 
modes (tangential to the boundary). Further polariza¬ 
tion effects arise with non-uniform metal coatings, since 
they have boundaries that are both tangential and radi¬ 
ally oriented relative to the cylindrical geometry of the 
fiber [T^ . It is therefore desirable to carry out two trans¬ 
mission measurements along the principal axis to observe 
directly such polarization effects on the strengths and 
positions of the cladding mode resonances [l9|. On the 
other hand, it has been shown here that a Jones ma¬ 
trix measurement can provide this information as well, 
in addition to other parameters of interest, such as PDF. 
While the PDF spectrum ’’hides” the physical effects re¬ 
sponsible for difference in transmission due to the differ¬ 
ent polarization states, it has been shown in the past to 
yield excellent limits of detection for surface plasmon res¬ 
onance based TFBG sensors We now proceed to com¬ 
pare the signal noise for refractive index measurements by 
the various polarization dependent data extraction tech¬ 
niques. 

As shown in Fig. 0 for a TFBG coated with a sparse 
layer of gold nanorods and immersed in water, the PDF 
parameter provides the absolute value of the difference 
between resonances observed in the transmission spectra 
measured along the principal axes. The relative position 
of the peaks, their amplitude, and their width become 
convoluted in the PDF parameter, which provides instead 
a maximum located somewhat in between the individual 
resonance maxima and a zero on either side correspond¬ 
ing to the wavelengths where the spectra cross each other. 


When the refractive index of the medium surrounding 


Pi .P2 , (dB) 




Wavelength, X (nm) 


FIG. 7. Two eigenvalue spectra (individual transmission 
along the principal axes) and corresponding polarization de¬ 
pendent loss (PDF) parameter. The peak in PDF spectrum 
is denoted by ”1” and zeros by ”2”. Here the more common 
dB scale is used. 

such TFBG changes the waveguiding characteristics of 
the cladding are modified and the resonances observed 
in the transmission spectrum change accordingly. There¬ 
fore, to detect changes in refractive index we can either 
follow the amplitudes and positions of individual reso¬ 
nances jl^ or of the PDF features. Here, the sensor 
was immersed in water and the refractive index was in¬ 
crementally increased in steps of An = 1.517 x 10-4 by 
adding 10 /i/ of ethylene glycol ( (727^4(OiJ )2 ) to 5 m/ 
to the water. The impact of each increase in refractive 
index on all parameters of interest is shown in Fi^IS] for 
a typical slice of the spectrum (it was noted in Q that 
there was little difference in wavelength shifts across the 
TFBG spectrum for this device). 

The refractive index change was chosen to have a rela¬ 
tively small value of (An = 1.517 x 10-4)to test the sen¬ 
sor detection limits and a linear fitting was used because 
the sensor response for such small changes is expected to 
be linear. The standard deviation of the errors from the 
linear fit was calculated in the usual manner by: 


SD = 


\ 


N 






(19) 


here /i is the expected value of x, and x = yappr — Vdata 
is the difference between the measured data ydata and its 
linear approximation yappr- 
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FIG. 8. Position of the resonance a) in the transmission spec¬ 
tra Ix and ly of light polarized along x and y geometrical axes 
here aligned with the principal axes, and b) in PDF spectra 
(the maximum and zero values of PDF are detected), as a 
function of refractive index change. The continuous lines rep¬ 
resent the least square approximation to the measured data, 
and SD is the standard deviation from the linear approxima¬ 
tion 

The results of the fits show that the most accurate 
detection of the refractive index change is achieved with 
the zeros in the PDF spectrum, as this measurement pro¬ 
vides the smallest standard deviation of SD = 1.40 pm. 
The worst result were obtained with the PDF peak 
{SD = 6.44 pm), while the detection of individual po¬ 
larized resonances provides an intermediate value of the 
standard deviation (essentially equal to 3 pm). It is not 
surprising that zeros of PDF should provide the most 
accurate results as they consist essentially of a differen¬ 
tial measurement between two spectra with well-defined 
crossing points that occur on the sides of the individ¬ 
ual resonances, where the spectral slope is highest. On 
the other hand, at the expense of an increase in noise 
by a factor of approximately 2, other effects such as dif¬ 
ferences in the change in the resonance amplitude for 


the two polarizations (which can be linked to differential 
loss or scattering) can be studied when the principal axis 
spectra are used. 


VI. CONCLUSION 

In this paper we investigated the use of Jones ma¬ 
trix and Stokes vector based techniques and polarization 
analysis to extract information from optical fiber sensors 
in non-polarization maintaining fibers. In particular we 
showed how to calculate transmission spectra with elec¬ 
tric fields E aligned with the system principal axes or 
along any other system axes. In the case of a TFBG 
inscribed in a non-polarization maintaining fiber for in¬ 
stance, the principal axes of the system are determined 
by the direction of the tilt of the grating planes (and its 
perpendicular). The transmission spectra for light po¬ 
larized in the tilt plane (P-polarized) or out of the tilt 
plane (S-polarized) can thus be extracted without hav¬ 
ing to separately align a linear polarizer upstream from 
the grating and hoping that the polarization remains lin¬ 
ear between the polarizer and the TFBG. Polarization- 
resolved spectra can also be obtained at a faster rate, 
for applications in chemical deposition process monitor¬ 
ing for instance, without the need to line up or rotate a 
polarizer between each measurement. 

The transmission spectra of P- and S-polarized light 
were compared with the TFBG spectral response along 
the principal axes. A small oscillation of about 8 degrees 
in the orientation of the principal axes as a function of 
wavelength was observed. 

Furthermore, it was determined that the best TFBG 
sensor detection limits are achieved when zeros in the 
PDF spectrum are followed, mainly because of the sharp¬ 
ness of crossing points occurring on the steep sides of the 
individual resonances, but at the expense of the direct ob¬ 
servation of the sensor response to polarized light (which 
are also available from the Jones matrix data). In the 
latter case, while the standard deviation of the spectral 
sensitivity is doubled relative to PDF measurements, ad¬ 
ditional information becomes available regarding the in¬ 
fluence of the measured medium on cladding mode loss, 
allowing further uses from the TFBG data sets [^ . 


[1] M. Zhang, K. Ihida-Stansbury, J. V. der Spiegel, and 
N. Engheta, “Polarization-based non-staining cell detec¬ 
tion,” Opt. Express 20 , 25378-25390 (2012). 

[2] G. A. Granch, G. M. H. Elockhart, and G. K. Kirk- 
endall, “Polarization properties of interferometrically in¬ 
terrogated fiber bragg grating and tandem-interferometer 
strain sensors,” J. Fightwave Technol. 24, 1787 (2006). 

[3] O. Kotov, M. Kol’chenko, and Y. E. Fozovik, “Ultra- 
high refractive index sensitivity of te-polarized electro¬ 
magnetic waves in graphene at the interface between two 


dielectric media,” Opt. Express 21 , 13533-13546 (2013). 

[4] R. Wang, F. Du, G. Zhang, Z. Man, Y. Wang, S. Wei, 
G. Min, S. Zhu, and X.-G. Yuan, “Plasmonic petal¬ 
shaped beam for microscopic phase-sensitive spr biosen¬ 
sor with ultrahigh sensitivity,” Opt. Fett. 38, 4770-4773 
(2013). 

[5] C. Gaucheteur, Y. Shevchenko, F.-Y. Shao, M. Wuil- 
part, and J. Albert, “High resolution interrogation of 
tilted hber grating spr sensors from polarization proper¬ 
ties measurement,” Opt. Express 19, 1656-1664 (2011). 











8 


[6] C. Caucheteur, S. Bette, C. Chen, M. Wuilpart, 
P. Megret, and J. Albert, “Tilted fiber bragg grating re- 
fractometer using polarization-dependent loss measure¬ 
ment,” Photonics Technology Letters, IEEE 20 , 2153- 
2155 (2008). 

[7] A. Bialiayeu, A. Bottomley, D. Prezgot, A. lanoul, and 
J. Albert, “Plasmon-enhanced refractometry using silver 
nanowire coatings on tilted fibre bragg gratings,” Nan¬ 
otechnology 23 , 444012-444022 (2012). 

[8] M. Z. Alam and J. Albert, “Selective excitation of ra¬ 
dially and azimuthally polarized optical fiber cladding 
modes,” J. Lightwave Technol. 31 , 3167-3175 (2013). 

[9] K. S. Lee and T. Erdogan, “Eiber mode coupling in trans¬ 
missive and reflective tilted fiber gratings,” Appl. Opt. 
39 , 1394-1404 (2000). 

[10] G. E. Villanueva, M. B. Jakubinek, B. Simard, C. J. 
Oton, J. Matres, L.-Y. Shao, P. Perez-Millan, and J. Al¬ 
bert, “Linear and nonlinear optical properties of car¬ 
bon nanotube-coated single-mode optical fiber gratings,” 
Opt. Lett. 36 , 2104-2106 (2011). 

[11] L.-Y. Shao, M. B. Jakubinek, T. Sun, B. Simard, and 
J. Albert, “Pour-wave mixing in carbon nanotube-coated 
optical fiber gratings,” Applied Physics Letters 100 , 
071108 (2012). 

[12] A. Bialiayeu, C. Caucheteur, N. Ahamad, A. lanoul, and 
J. Albert, “Self-optimized metal coatings for fiber plas- 
monics by electroless deposition,” Optics Express 19 , 
18742-18753 (2011). 


[13] R. C. Jones, “A new calculus for the treatment of optical 
systems,” J. Opt. Soc. Am. 31 , 488-493 (1941). 

[14] B. L. Heffner, “Deterministic, analytically complete 
measurement of polarization-dependent transmission 
through optical devices,” Photonics Technology Letters, 
IEEE 4, 451-454 (1992). 

[15] H. G. Berry, G. Gabrielse, and A. E. Livingston, “Mea¬ 
surement of the stokes parameters of light,” Appl. Opt. 
16 , 3200-3205 (1977). 

[16] D. H. Goldstein, Polarized Light, Revised and Expanded 
(CRC Press, 1992). 

[17] W. H. McMaster, “Polarization and the stokes parame¬ 
ters,” American Journal of Physics 22 , 351-362 (1954). 

[18] W. Zhou, D. J. Mandia, M. B. Griffiths, A. Biali¬ 
ayeu, Y. Zhang, P. G. Gordon, S. T. Barry, and J. Al¬ 
bert, “Polarization-dependent properties of the cladding 
modes of a single mode fiber covered with gold nanopar¬ 
ticles,” Opt. Express 21 , 245-255 (2013). 

[19] C. Caucheteur, C. Chen, V. Voisin, P. Berini, and J. Al¬ 
bert, “A thin metal sheath lifts the eh to he degeneracy 
in the cladding mode refractometric sensitivity of opti¬ 
cal fiber sensors,” Applied Physics Letters 99 , 041118 
( 2011 ). 

[20] W. Zhou, D. J. Mandia, M. B. E. Griffiths, S. T. Barry, 
and J. Albert, “Effective permittivity of ultrathin chemi¬ 
cal vapor deposited gold films on optical fibers at infrared 
wavelengths,” The Journal of Physical Chemistry C 118 , 
670-678 (2014). 



